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Abstract 

First, we show how to reduce the sensitivity of the NLO predictions of the Drell- 
Yan production of low-mass, lepton-pairs, at high rapidity, to the choice of factorization 
scale. In this way, observations of this process at the LHC can make direct measurements 
of parton distribution functions in the low x domain; x < 10~ 4 . Second, we find an 
inconsistency in the conventional NLO treatment of the infrared region. We illustrate the 
problem using the NLO coefficient function of Drell-Yan production. 



1 LHC Drell-Yan production as a probe of low x 

The very high energy of the LHC allows a probe of the parton distribution functions (PDFs) 
of the proton at extremely small x, a region not accessible at previous accelerators. To extract 
the PDFs we describe the experimentally observed cross sections as a convolution of the PDFs 
and the cross section for the hard partonic subprocess, which is of the form 

da/d 3 p = J dx x dx 2 PDF(x 1 , fi F ) \M(p;fi F ,fi R )\ 2 PDF(x 2 ,/i F ) , (1) 

where a sum over the various pairs of PDFs is implied. 

Here we focus on Drell-Yan production of a low mass fi + fi~ pair. At LO the production 
of a ji + ji~ system of mass M and rapidity Y arises from the subprocess 7* — > qq with x lj2 = 



M exp(±Y) / y/s . So for M = 6 GeV, Y = 4 (a domain which is accessible to the LHCb 
experiment) we probe x\ = 4.7 x 1CT 2 , xi = 1.6 x 1CT 5 . The problem is that, in the low x 
region, the PDFs strongly depend on the choice of the factorization scale fi F , see Fig. [l|a). It 
is made worse due to the dominance of the gluon PDF at small x, which means that the LO 
qq — > 7* subprocess is overshadowed by the NLO subprocess gq — > q-y*- However, it is this 
very dominance which will allow us to introduce a procedure which greatly suppresses the scale 
dependence of the predictions. 
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Figure 1: (a) Sensitivity of M = 6 GeV Drell-Yan production, as a function of rapidity Y, 

to the choice of factorization scale: \i F = M/2, M, 2M, at LO, NLO, NNLO. (b) The bold lines 
correspond to the choice fi F = /i = 1AM which minimizes C^°, and show the stability with 
respect to the variations pL = M/2, M, 2M in the scale of the PDFs convoluted with - the 

li dependence is indicated by the symbolic equation at the top of the diagram. The dashed lines 
show that the stability disappears for other choices of fi . The small crosses are the NNLO result. 
In this figure, taken from [lj, the renormalization scale is fixed at = M. 



The plan is to choose the value of \xp which minimizes the higher order a s NLO, NNLO,.. 
contributions. To sketch the idea, we start with the LO expression for the cross section: 
cr(/ip) = PDF(/ip) <g> C h ° ® PDF(/ip) in the collinear approach. The effect of varying the scale 
from m to fip, in both the left and right PDFs, can be expressed, to first order in a s , as 



a(p F ) = PDF(m) ® (C LO + gin ( (P lcft C LO + C LO P rig ht)J <8> PDF(m), (2) 

where the splitting functions P r j g ht = P qq + P qg and P\ e & = Pqq + Pq g act on the right and left 
PDFs respectively. We may equally well have incoming <fs in P r i g ht and incoming g's in Pi e ft. 

At NLO we may write 

a(fi F ) = PDF(fi F )® (C L0 + a s C c N o ^°)®PDF(/i F ), (3) 
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where the 2 — > 2 subprocesses qq — > g*y* and gq — )■ g7* are now calculated with better, than 
LLA, accuracy. However, we must subtract from C NLO , the part of the contribution already 
included, to LLA accuracy, in the a s term in The remaining contribution, (^^(/i^), now 
depends on the scale /ip, coming from the ^f dependence of the LO LLA term that has been 
subtracted off. The trick is to choose an appropriate scale, fip = fiQ, so as to minimize the 
remaining NLO contribution (/Af)- Explicit calculation [TJ shows that the optimum choice 
is \ip — A*o = 1-4M. The stability of the prediction, using MSTW NLO PDFs [2], is shown in 
Fig. [ijb). For Y > 3, pure DGLAP PDF extrapolations become unreliable due to the absence 
of absorptive, ha(l/x),.. modifications. Rather, LHC data will provide a direct measure of PDFs 
in this low x domain. 



2 Treatment of the infrared region in perturbative QCD 

Interestingly, a spin-off of the above study highlighted an inconsistency in the conventional 
treatment of the infrared region [3]. Again we use Drell-Yan as an example. For the main NLO 
subprocess we have 



— > g7*) a 2 a s z 1 



d\t\ 9M 2 \t\ 



(4) 



where z = M 2 /s and Vs is the incoming gq cm. energy. (Strictly speaking, z is the ratio 
of the light-cone momentum fraction carried by the 'daughter' quark to that carried by the 
'parent' gluon, z = x^/x^.) In order to calculate the inclusive cross section da/dM 2 , it seems 
that we have to integrate over t starting from t = 0. If this were necessary, then we would face 
an infrared divergency. 

However, we follow the procedure of the previous Section, which we call the physical ap- 
proach. To avoid double-counting, we subtract the LO DGLAP-generated contribution. Then 
the remaining contribution of Q is [3] 



ddfj^igq— ^97*) a 2 a s z 1 |\^ ,o o, o t 2 _ „ t 



d\t\ ~ 9M 2 \t\ 



[{1-zY + z 2 } 0(1*1-/4) + z 2 — - 2z- Af2 



(5) 



which has no singularity as t — > 0. The LO DGLAP evolution has accounted for all virtualities 
\t\ = k 2 < [i 2 F \ with the contribution of k 2 < Qq hidden in the phenomenological input PDF at 

Q 2 = Ql 

On the other hand, the conventional prescription evaluates the inclusive cross section, 
da/dM 2 , by integrating Q over the infrared divergency at t = using 4 + 2e dimensional 
space to regularize the integral [31 El E] - Then the contribution from very small t produces a 
1/e pole, which is absorbed into the incoming PDF. The conventional prescription is as follows. 
The same gq — > q^* diagram, but now generated by LO DGLAP evolution, is considered; it 
gives an 1/e pole which cancels the corresponding 1/e pole in hard matrix element (coefficient 
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function). However, we are left with e/e terms of infrared origin, which produce a non-zero 
result as e — > 0. Unlike the finite e/e terms of ultraviolet origin, which can be treated as point- 
like counter-terms in the Lagrangian, the infrared e/e contribution makes no physical sense in 
QCD theory, since the confinement eliminates any interaction at very large distances. 

To be explicit, the conventional prescription gives 

M 2 da^(gq^ qi *) a 2 a s z IV , 2 , 21 , (1 - zf , 1 , Q 7 



|[(l-.) 2 + . 2 ]lnil-^ + i + 3,-^ 2 }, (6) 



dM 2 9M 2 

which is to be compared with the result, 

M'da^jgq -> qj*) _ a^z ( { _ 2 2] (1 - z) 1 _ 3 ,1 

ciM 2 ~ 9M 2 \ Ll Z J +*J m 3 + 2 +2; 2 ^J> 10 

obtained, from ([5]), using the physical approach. We have checked that the spurious e/e contri- 
bution in the conventional approach is responsible for the difference, 

^{[(l-.) 2 + . 2 ]ln(l-.) + 2,(l-.)}, (8) 

between ^ and ([7]). 

The message is that to get a reliable result we must use the physical approach and to subtract 
from ^ the 1/t singularity exactly, so obtaining the non-singular expression ([5]). Dimensional 
regularisation is not appropriate in the infrared region. As an example we have used the NLO 
coefficient function for Drell-Yan production. However, the result applies more generally. It 
is relevant to other processes, such as deep inelastic scattering [3]. It is relevant to the NLO 
splitting functions. 

It should be noted that the physical approach cannot be considered as an alternative fac- 
torization scheme. That is, the difference between the conventional and the physical coefficient 
functions cannot be compensated by a re-definition of the parton distributions (PDFs). The 
corresponding (infrared) e/e corrections in the splitting functions do not coincide with those 
that come from the re-definition of the PDFs. 
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